" esp@cetiet document view Page 1 of 1 



Method for determining density distributions and atomic number 
distributions during radiographic examination methods 



Also published as: 

'1 WO03024331 (A3) 
'1 WO03024331 (A2) 
'1 US2004223585 (A1) 

Classification: 

- international: Q01 N23/083; A61 B6/00 

- european: A61 B6/00D; G01 N23/04 
Application number: DE2001 1 0431 31 20010903 
Priority number(s): DE2001 1 0431 31 2001 0903 



Patent number: DE1 01431 31 
Publication date: 2003-04-03 

inventor: HEISMANN BJOERN (DE); STIERSTORFER KARL 

(DE) 

Applicant: SIEMENS AG (DE) 



Report a data error here 



Abstract of DE10143131 

The invention relates to a method, an X-ray apparatus and a computer software product, which make it 
possible to determine the distributions of material density and atomic number in objects to be examined of 
the human or animal type and in objects to be examined that are used in material and safety engineering. 
The functional dependency of a first X-ray absorption value of density and atomic number is determined in 
the instance of a first X-ray spectrum, and at least the functional dependency of a second X-ray 
absorption value of density and atomic number is determined in the instance of a second X-ray spectrum. 
Based on a recording of a first distribution of X-ray absorption values of the object to be examined in the 
instance of a first X-ray spectrum, and on a recording of at least one second distribution of X-ray 
absorption values of the object to be examined in the instance of a second X-ray spectrum, the values for 
density and atomic number are determined by comparing the functional dependency of a first X-ray 
absorption value of the first distribution of X-ray absorption values with the functional dependency(ies) of 
the X-ray absorption values, which are assigned to the first X-ray absorption value, of the second and/or 
other distributions of X-ray absorption values. 
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Prufungsantrag gem. § 44 PatG ist gestellt 

® Verfahren zur Ermittlung von Dichte- und Ordnungszahlverteilungen bei radiographischen 
Untersuchungsverfahren 

(57) Beschrieben werden ein Verfahren, eine Rontgenappa- r 
ratur und ein Computersoftwareprodukt, die ein Ermitteln 
der Verteilungen von Materialdichte und Ordnungszahl in 
Untersuchungsobjekten menschlicher wie tierischer Na- 
tur und in Untersuchungsobjekten der Material- und Si- 
cherheitstechnik ermoglichen. Es wird die funktionale Ab- 
hangigkeit eines ersten Rontgenabsorptionswertes von 
Dichte und Ordnungszahl bei einem ersten Rontgen- 
strahlspektrum und zumindest die funktionale Abhangig- 
keit eines zweiten Rontgenabsorptionswertes von Dichte 
und Ordnungszahl bei einem zweiten Rontgenstrahlspek- 
trum bestimmt. Ausgehend von einem Aufzeichnen einer 
ersten Verteilung von Rontgenabsorptionswerten des Un- 
tersuchungsobjekts bei einem ersten Rontgenstrahlspek- 
» trum und einem Aufzeichnen zumindest einer zweiten 
i Verteilung von Rontgenabsorptionswerten des Untersu- 
chungsobjekts bei einem zweiten Rontgenstrahlspektrum 
werden die Werte fur Dichte und Ordnungszahl aus einem 
Vergleich der funktionalen Abhangigkeit eines ersten 
Rontgenabsorptionswertes der ersten Verteilung von 
Rontgenabsorptionswerten mit der oder den funktionalen 
Abhangigkeiten der dem ersten Rontgenabsorptionswert 
zugeordneten Rontgenabsorptionswerte der zweiten und/ 
oder weiteren Verteilungen von Rontgenabsorptionswer- 
ten ermittelt. 




BUNDESDRUCKEREI 02.03 103 140/259/1 
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Beschreibung 

[0001] Die voriiegende Erfindung betrifft ein Verfahren, eine Rontgenapparatur und ein Computersoftwareprodukt zur 
Ermittlung von Bilddaten tfber die Verteilung physikalisch und chemisch relevanter Daten in Untersuchungsobjekten 
5 menschlicher oder tierischer Natur oder aus dem Bereich der Material- oder Sicherheitspriifung. Insbesondere betrifft die 
Erfindung ein Auftrennen radiographischer Bilddaten in eine Verteilung der Materialdichte und eine Verteilung der Ord- 
nungszahl im Untersuchungsobjekt 

[0002] Ergebnis aller radiographischen Verfahren, wie beispielsweise der Cornputer-Tomographie, der Mammogra- 
phie, der Angiographic, der Rontgen-Inspektionstechnik oder vergleichbarer Verfahren, ist die Darstellung der Schwa- 
10 chung eines Rontgenstrahls entlang seines Weges von der Rontgenquelle zum Rontgendetektor. Diese Sch wachung wird 
von den durchstrahlten Medien bzw. Materi alien entlang des Strahlengangs verursacht. Die Sch wachung wird Ublicher- 
weise als der Logarithmus des Verhaltnisses der Intensitat der geschwachten zur Primarstrahlung definiert und bezogen 
auf ein Wegnormal als Schwachungskoeffizient des Materials bezeichnet 

[0003] Eine Vielzahl radiographischer Untersuchungsgerate verwendet zur Darstellung der Schwachungsverteilung ei- 
15 nes Rdntgenstrahls in einem Untersuchungsgegenstand statt des Schwachungskoeffizienten einen auf den Schwachungs- 
koefflzienten von Wasser normierten Wert, die CT-Zahl. Diese berechnet sich aus einem aktuell durch Messung ermittel- 
ten Schwachungskoeffizienten u und dem Referenz-Schwachungskoeffizienten u H2 o nach folgender Gleichung: 



20 



35 



C = 1000 X [HU] (1) 

Mh,o 



mit der CT-Zahl C in der Einheit Hounsfield [HU]. Fur Wasser ergibt sich ein Wert Ch 2 o = 0 HU und fur 0 Luft ein Wert 
C L = -1000 HU. 

[0004] Da beide Darstellungen ineinander transformierbar bzw. Equivalent sind, bezeichnet im folgenden der allge- 
25 mein gewahlte Begriff Schwachungswert sowohl den Schwachungskoeffizienten u als auch den CT-Wert. Femer werden 
im Sachzusammenhang dieser Erfindungsbeschreibung die Begriffe Material und Gewebe austauschbar verwendet. Es 
wird unterstellt, dass ein Material im Kontext einer medizinisch angezeigten Untersuchung ein anatomisches Gewebe 
sein kann, und urngekehrt unter Gewebe in der Material- und Sicherheitspriifung ein beliebiges Material eines Untersu- 
chungsobjekts zu verstehen isL 
30 [0005] Erhohte Schwachungswerte lassen sich entweder auf Materialien hoherer Ordnungszahl, wie beispielsweise 
Calcium im Skelett oder Jod in einem Kontrastmittel zuriickfuhren, oder auf eine erhohte Weichteildichte, wie etwa bei 
einem Lungenknoten. Der lokale Schwachungskoeffizient u am Ort 7 ist abhangig von der in das Gewebe bzw. Material 
eingestrahlten Rontgenenergie E und der iokalen Gewebe- bzw. Materialdichte p entsprechend der folgenden Gleichung: 



■(5) 



d = M(E,r) = - E)xp (r) (2) 



mit dem energie- und materialabhangigen Massenschwachungskoeffizienten 




45 [0006] Die energieabhangige Rontgenabsorption eines Materials, wie sie von seiner effektiven Ordnungszahl be- 
stimmt wird, iiberlagert daher die von der Materialdichte beeinflusste Rontgenabsorption. Materialien bzw. Gewebe un- 
terschiedlicher chemischer wie physikalischer Zusammensetzung konnen daher im Rontgenbild identische Schwa- 
chungswerte aufweisen. Urngekehrt kann dagegen aus dem Schwachungswert einer Rontgenaufhahme nicht auf die Ma- 
terialzusammensetzurig eines Untersuchungsobjekts geschlossen werden. 

50 [0007] Eine korrekte Interpretation der somit eigentlich eher unanschaulichen Verteilung der Schwachungswerte in ei- 
nem mit einem radiographischen Untersuchungsverfahren erstellten Rontgenbild kann im medizinischen Sektor meist 
nur aufgrund morphologischer Kriterien erfolgen und erfordert meist einen Radiologen mit jahrzehntelanger Erfahrung 
auf seinem Gebiet Dennoch konnen in einigen Fallen Strukturen, die in der Bildgebung einer Rontgenuntersuchung mit 
erhohten Schwachungswerten auffallen, nicht klar klassifiziert werden. Beispielsweise ist eine hilusnahe Verkalkung auf 

55 einer Thoraxiibersichtsaufhahme nur schwer von einem orthograd zur Bildebene liegenden GefaB zu unterscheiden. 
Auch eine diffuse Kalkeinlagerung kann beispielsweise kaum von einer frischen Einblutung unterschieden werden. 
[0008] Auch in der Material- und Sicherheitspriifung erganzt der Prufer i. a. die Information der Darstellung einer 
Schwachungswert- Verteilung durch seine personliche Fachkenntnis und berufliche Erfahrung. Dennoch ist ihm z. B. ein 
sicheres Unterscheiden einer kunststoffgebundenen Sprengstoffmischungen von einem nichtexplosiven Kunststoff di- 

60 rekt aus einem Rontgenbild nicht moglich. 

[0009] Hierfur sind Verfahren zur Darstellung materialcharakteristischer Werte erforderlich^^^^alej^aget. al be- 
schreiben in "Materialselektive Bildgebung und Dichtemessung mit der Zwei-Spektren-Methode, L Grundlagen und 
Methodik, W. Kalender, W. Bautz, D. Felsenberg, C. SuB und E. Klotz, Digit. Bilddiagn. 7, 1987, 66-77, Georg Thieme 
Verlag" ein Verfahren zur Basismaterialzerlegung bei Rontgenaufhahmen. Das Verfahren basiert auf dem EfFekt, dass 

65 Materialien und Gewebe hoherer Ordnungszahl niederenergetische Rontgenstrahlung deutlich starker absorbieren als 
Materialien bzw. Gewebe niederer Ordnungszahl. Bei hoheren Rontgenstrahlenergien gleichen sich dagegen die Schwa- 
chungswerte an und sind vorwiegend eine Funktion der Materialdichte. 

[0010] Im Kontext dieser Beschreibung wird der Begriff Ordnungszahl, soweit nicht anders angegeben, nicht im stren- 
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gen, elementbezogenen Sinn verwendet, sondem bezeichnet stattdessen eine effektive Ordnungszahl eines Gewebes, re- 
spektive Materials, die sich aus den chemischen Ordnungszahlen und Atomgewichten der am Aufbau des Gewebes bzw. 
Materials beteiligten Elemente berechnet 

[0011] Im von W. Kalender et. al vorgeschlagenen Verfahren werden die Rontgenschwachungswerte eines Untersu- 
chungsobjekts mit Rontgenstrahlen niederer und hoherer Eneigie gemessen und die erhaltenen Werte mit den entspre- 5 
chenden Referenzwerten zweier Basismaterialien wie beispielsweise Calcium (fur Skelettmaterial) und Wasser (fur 
Weichteilgewebe) verglichen. Es wird angenommen, dass sich jeder Messwert als lineare Superposition der Messwerte 
der beiden Basismaterialien darstellen lasst Zum Beispiei kann fur jedes Element der bildlichen Darstellung des Unter- 
suchungsobjekts aus dem Vergleicb mit den Werten der Basismaterialien ein Skelettanteil und ein Weichgewebeanteil 
berechnet werden, so dass eine Transformation der ursprunglichen Aufhahmen in Darstellungen der beiden Basismate- to 
rialien Skelettmaterial und Weichteilgewebe resultiert. 

[0012] Die Basismaterialzerlegung bzw. das Zwei-Spektren- Verfahren eignet sich damit zur Auftrennung bzw. Unter- 
scheidung von anatomischen Strukturen in menschlichen und tierischen Geweben mit stark unterschiedlicher Ordnungs- 
zahl. In der Material- und Sicherheitspriifung konnte damit beispielsweise eine Auftrennung nach vordefinierten Mate- 
rialarten, sogenannten Materialklassen, erfolgen. Eine funktionelle Darstellung, die physikalische und chemische Cha- 15 
rakteristiken der untersuchten Materialien bzw. Variationen dieser Charakteristiken innerhalb einer Materialart erkennen 
lasst, ist nicht Zielrichtung der Basismaterialzerlegung. 

[0013] Es ist daher Aufgabe der vorliegenden Erfindung ein Verfahren, eine Rontgenapparatur und ein Computersoft- 
wareprodukt vorzuschlagen, die auf der Basis radiographisch erstellter Absorptionsdaten eines Untersuchungsobjekts 
chemisch und physikalisch relevante Information iiber das Untersuchungsobjekt erstellen. 20 
[0014] Diese Aufgabe wird gelost durch ein Verfahren zur Ermittlung der Verteilungen von Dichte und Ordnungszahl 
in einem Untersuchungsobjekt mit Schritten zum Aufzeichnen einer ersten Verteilung von Rontgenabsorptionswerten 
des Untersuchungsobjekts bei einem ersten Rontgenspektrum und zum Aufzeichnen zumindest einer zweiten Verteilung 
von Rontgenabsorptionswerten des Untersuchungsobjekts bei einem zweiten Rontgenstrahlspektrum, zum Bestimmen 
einer ersten funktionalen Abhangigkeit eines ersten Rontgenabsorptionswertes der ersten Verteilung von Rontgenab- 25 
sorptionswerten von Dichte und Ordnungszahl und zum Bestimmen einer zweiten funktionalen Abhangigkeit eines, dem 
ersten Rontgenabsorptionswert zugeordneten zweiten Rontgenabsorptionswertes der zweiten Verteilung von Rontgenab- 
sorptionswerten von Dichte und Ordnungszahl, und zum Ermitteln eines Werts fur die Dichte und die Ordnungszahl aus 
einem Vergleich der ersten funktionalen Abhangigkeit mit der zweiten funktionalen Abhangigkeit und/oder weiteren 
funktionalen Abhangigkeiten. 30 
[0015] Der im Zusammenhang dieser Schrift verwendete Begriff "Rontgenspektrum" besitzt eine weiter gefasste Be- 
deutung als nur die Spektralverteilung einer von der Rontgenquelle der Apparatur emittierten Rontgenstrahlung. Auf 
Seiten der Rontgendetektoren werden unterschiedliche Spektralanteile einer Strahlung mit unterschiedlichen Wirkungs- 
graden umgesetzt und somit verschieden gewichtet. Die daraus resultierende effektive Spektralverteilung wird in dieser 
Schrift als Rontgenspektrum bezeichnet. 35 
[0016] Obige Aufgabe wird weiterhin gelost durch eine Rontgenapparatur zur Ermittlung der Verteilungen von Dichte 
und Ordnungszahl in einem Untersuchungsobjekt mit einer Rontgenquelle zur Emission von Rontgenstrahlung, einem 
Rontgendetektor zum Nachweis der von der Rontgenquelle emittierten Rontgenstrahlung und zum Umsetzen der Ront- 
genstrahlung in elektrische Signale fur eine weitere Verarbeitung, und einer Signalverarbeitungseinrichtung fur das Ver- 
arbeiten der elektrischen Signale des Rontgendetektors, wobei die Signalverarbeitungseinrichtung die Verteilungen von 40 
Dichte und Ordnungszahl in einem Untersuchungsobjekt auf Basis zumindest zweier, an der Rontgenapparatur mit un- 
terschiedlichen Rontgenstrahlspektren aufgezeichneten Verteilungen einer Rontgenabsorption des Untersuchungsob- 
jekts nach einem erfindungsgemaBen Verfahren ermittelt. 

[0017] Femer wird die Aufgabe gelost durch ein Computersoflwareprodukt zur Ermittlung der Verteilungen von 
Dichte und Ordnungszahl in einem Untersuchungsobjekt auf Basis von zumindest zwei bei unterschiedlichen Rontgen- 45 
strahlspektren aufgezeichneten Bilddatensatzen der Verteilung von Rontgenabsorptionswerten in dem Untersuchungsob- 
jekt nach einem erfindungsgemaBen Verfahren. 

[0018] Die vorliegende Erfindung ermoglicht die Berechnung der raumliche Verteilung der mittleren Dichte p(7) und 
der effektiven Ordnungszahl Z(7) aus einer Auswertung der spektral beeinflussten Messdaten einer Rontgenapparatur. 
Man erhalt hieriiber neuartige Kontraste, insbesondere beziiglich der chemischen und physikalischen Zusammensetzung 50 
des Untersuchungsobjekts. Diese, bislang Magnetresonanz-Systemen vorbehaltene funktionale Darstellung eines Unter- 
suchungsobjekts eroffhet der Rontgendiagnostik wie der Rontgeninspektionstechnik eine Vielzahl neuer Anwendungen. 
[0019] Beispielsweise eriaubt die Darstellung der Verteilung der Ordnungszahl im Gewebe u. a. Einblicke in die bio- 
chemische Zusammensetzung eines untersuchten Objekts, Kontraste aufgrund des chemischen Aufbaus in bisher dichte- 
homogen dargestellten Organen, eine quantitative Bestimmung von Korperbestandteilen wie z. B. Jod oder dergleichen 55 
und ein Heraussegmentieren von Calcifizierungen basierend auf der Ordnungszahl. Die isolierte Dichtedarstellung eines 
Objekts gestattet eine genaue Schwerpunktbestimmung und Dichtevermessung von Objekten, wie sie u. a. beispiels- 
weise bei Osteoporose vorgenommen werden. 

[0020] Im Bereich der Sicherheitstechnik bedeutet dies eine zuverlassigere Detektierbarkeit gefahrlicher Komponen- 
ten, insbesondere von Explosivstoffen. In der Materialpriifung eroffnet sich der Zugang zur quantitativen Untersuchung 60 
der Materialzusammensetzung und der Dichteverteilung in den Priiflingen. 

[0021] Weitere vorteiihafte Ausfiihrungen der vorliegenden Erfindung sind in den entsprechenden Unteranspriichen 
definiert. 

[0022] Vorteilhaft erfolgt die Bestimmung der funktionalen Abhangigkeit der Rontgenabsorptionswerte von Dichte 
und Ordnungszahl fur zumindest ein Rontgenstrahlspektrum mittels Referenzmessung an einer Eichprobe und/oder in 65 
Form einer Simulation auf der Basis eines physikalischen Modells, so dass eine anlagenspezifische Abhangigkeit eines 
Rontgenabsorptionswertes von Dichte und Ordnungszahl erhalten wird. 

[0023] Vorzugs weise wird ein Umformen der Verteilungen der Rontgenabsorptionswerte in eine Verteilung der Dichte 
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und eine Verteilung der Ordnungszahl ftir jeden der zugeordneten Rontgenabsorptionswerte der ersten \ferteilung der 
R6ntgenabsorpdonswerte und der weiteren Verteilungen der Rontgenabsorptionswerte auf der Grundlage der Ermittlimg 
eines Wertepaares ftir Dichte und Ordnungszahl so vorgenommen, dass das Wertepaar die bestirnmten funktionalen Ab- 
hangigkeiten der Rontgenabsorption von Dichte und Ordnungszahl fur das erste Rontgenstrahlspektrum und zumindest 
5 ein weiteres R6ntgenstrahlspektrum erfullt. Damit konnen Dichte und Ordnungszahl fur ein Bildelement einfach als 
Schnittmenge der funktionalen Abhangigkeiten der einander zugeordneten Rontgenabsorptionswerte der aufgezeichne- 
ten Verteilungen der Rontgenabsorptionswerte berechnet werden. 

[0024] Vorteilhafterweise weist das erste Rontgenspektrum eine Quantenenergie auf, die relativ zur Quantenenergie 
des zweiten Rontgenspektrums eine Rontgenabsorption durch den Photoeffekt begiinstigt, so dass eine hohe Auflosung 

to in der Bestimmung der Ordnungszahlen erhalten wird. 

[0025] In einer bevorzugten Ausfuhrungsform der vorliegenden Erfindung wird zurn Verandem eines Rontgenstrahl- 
spektrums fUr das Aufzeichnen des Untersuchungsobjekts eine Veranderung zumindest eines Betriebsparameters der 
Rdntgenrohre vorgenommen, wobei die Rontgenquelle in einem ersten Betriebszustand ein erstes Rontgenstrahlspek- 
trum und in einem zweiten Betriebszustand ein davon verschiedenes zweites Rontgenstrahlspektrum emittiert, so dass 

15 ein schneller Wechsel zwischen zwei Rontgenspektren ermbglicht ist. 

[0026] Ferner wird zum Verandem eines Rontgenstrahlspektrums fur das Aufzeichnen des Untersuchungsobjekts vor- 
teilhaft eine Veranderung der Detektorcharakteristik vorgenommen, wobei der R6ntgendetektor spektrale Teilbereiche 
der von der RSntgenqueile empfangenen Rontgenstrahlung in voneinander unabhangige elektrische Signale umsetzt und 
hierbei ein simultanes Aufzeichnen von Verteilungen der Rontgenabsorption bei unterschiedlichen Rontgenspektren zu- 

20 lasst, 

[0027] Das Aufzeichnen zumindest einer Teilregion des Untersuchungsobjekts erfolgt ferner vorteilhaft mit erhohter 
Rontgendosis damit fur diese Teilregion eine zufriedenstellende Auflosung der Dichte- und/oder Ordnungszahl-\fertei- 
lung erzielt wird ohne das gesamte Untersuchungsgebiet mit einer erhohten Rontgendosis zu belegen und z. B. einen Pa- 
tienten einer zu hohen Strahlenbelastung auszusetzen. 

25 [0028] Zum Aufzeichnen einer ersten Verteilung einer Rontgenabsorption des Untersuchungsobjekts kann entspre- 
chend einer vorteilhaften Weiterentwicklung ein Material in den Strahlengang der Rontgenstrahlung eingebracht werden 
das zum Aufzeichnen einer zweiten Verteilung einer Rontgenabsorption des Untersuchungsobjekts nicht in den Strahlen- 
gang der Rontgenstrahlung eingebracht wird. Damit konnen beim Aufzeichnen der ersten Schwachungswertverteilung 
vor bzw. nach dem Untersuchungsobjekt Spektralanteile der Rontgenstrahlung, die fur das jeweils verwendete Material 

30 spezifisch sind, herausgefiltert werden, so dass eine hohe Empfindlichkeit im Nachweis einer, diesem Material entspre- 
chenden Substanz im Untersuchungsobjekt erreicht wird. 

[0029] Im Folgenden wird die vorliegende Erfindung anhand von Ausfuhrungsbeispielen naher beschrieben, wobei auf 
folgende Figuren verwiesen wird, von denen 

[0030] Fig. 1 anhand einer Isoabsorptionslinie das Zustandekornmen identischer Schwachungswerte u bei Materialien 
35 unterschiedlicher Zusammensetzung zeigt, 

[0031] Fig. 2 die Energieabhangigkeit der Rontgenschwachung fur drei Elemente zeigt, 

[0032] Fig. 3 a ein Funktionsschema der erfindungsgemaBen Berechnungsmethoden zur Ermittlung von Isoabsorpti- 
onsUnien zeigt, 

[0033] Fig. 3b ein Ablaufdiagramm der Transformation der Rontgenabsorption werte in Werte der Materialdichte und 
40 Ordnungszahl zeigt, 

[0034] Fig. 4 zwei Isoabsorptionslinien einer Gewebeart bei zwei unterschiedlichen Rontgenstrahlspektren zeigt, 
[0035] Fig. 5 schematisch drei typische Rontgentechniken zeigt und 

[0036] Fig. 6 schematisch verschiedene Methoden zur Variation der Rontgenstrahlspektren darstellt. 

[0037] Die Isoabsorptionslinie 11 des Diagramms 10 der Fig. 1 verbindet alle Wertepaare (p, Z) mit, bei einem defi- 

45 nierten Rontgenstrahlspektrum identischem Schwachungswert u bzw. C. Die Darstellung der Fig. 1 verdeutlicht, dass In- 
formation iiber Art und Zusammensetzung eines Gewebes bzw. Materials nicht allein auf die Schwachungswerte eines 
Rontgenbildes gestiitzt abgeleitet werden konnen. Oblicherweise benutzt ein Radiologe zum Identifizieren von Gewebe- 
arten im Rontgenbild seine Anatomiekenntnisse und sucht auf dieser Basis nach UnregelmaBigkeiten. Zur Klarung der 
Identitat der UnregelmaBigkeiten ist ein Mediziner dann wiederum gezwungen auf Erfahrungs werte und morphologische 

50 Kriterien zuriickgreifen. Ahnlich stutzt sich ein Fachmann der Material- uns Sicherheitsprufung zur Beurteilung des ra- 
diographischen Befunds auf seinen beruflichen Erfahrungsschatz. 

[0038] Rontgenstrahlung wird von unterschiedlichen Materialien und abhangig von der Energie der Rontgenstrahlung 
unterschiedlich stark geschwacht. Die Fig. 2 illustriert dies anhand der Energieabhangigkeit 20 des Massenschwa- 
chungskoeffizienten fur Wasser 23, Calcium 22 und Jod 21. Dies ist auf unterschiedlich wirkende Schwachungsmecha- 

55 nismen bei den verschiedenen Materialien zuriickzufuhren. Im diagnostisch relevanten Energiebereich der Rontgen- 
strahlung ist die Rontgenschwachung im wesentlichen auf die Absorption, verursacht durch den Photoeffekt und die, auf 
dem Comptoneffekt beruhende Streuung zuriickzuftihren. Die Absorption ist insbesondere bei niedriger Energie der 
Rontgenstrahlung und bei Geweben mit hoher Ordnungszahl relevant. Die Streuung weist eine geringe Abhangigkeit 
von der Energie der Rontgenstrahlung auf und ist im wesentlichen von der Elektronendichte, vermittelt iiber die physi- 

60 kalischen Dichte des Gewebes, abhangig. 

[0039] Die im Kontext dieser Beschreibung vereinf acht als Ordnungszahl titulierte effektive Ordnungszahl Z einer be- 
stirnmten Gewebeart errechnet sich aus den Ordnungszahlen Zi der am Aufbau beteiligten Elemente, deren Atomge- 
wichte Aj und deren lokalen materialaquivalenten Dichten pi beispielsweise zu: 

65 
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[0040] FUr reines Calcium erhalt man Zca = 20, fur Calciumhydrid ca. Zcan 2 = 16,04 und fiir Wasser etwa Z H2 o = 
7,428. Die chemische oder auch biochemische Zusammensetzung eines Objekts kann daher sehr gut uber die Ordnungs- 
zahl Z erfasst werden. 

[0041] Voraussetzung fiir eine Berechnung der Ordnungszahl- und Dichteverteilung in einem Untersuchungsgebiet 
sind zumindest zwei, in der Aufnahmegeometrie identische, aber mit unterschiedlicher Energie der angewandten Ront- 15 
genstrahlung erstellte Rontgenaufnahmen des Gebiets. Bei Verwendung von mehr als zwei mit unterschiedlicher Ront- 
genstrahlenergie aufgezeichneten Rontgenaufnahmen konnen die Z- und p-Auflosung verbessert werden, doch erhoht 
sich dadurch auch die Strahlenbelastung. Im Falle der Untersuchung eines Patienten ist diese Moglichkeit daher nicht im- 
mer gegeben. 

[0042] Ausgangspunkt der Umwandlung von Schwachungswert basierenden Bilddaten in Verteilungsbilder der Ord- 20 
nungszahlen und der Material- bzw. Gewebedichte ist die Kenntnis der Isoabsorptionslinien fur jedes Rontgenspektrum 
einer Rontgenapparatur. Wie bereits erwahnt, ist hierbei unter Rontgenspektrum nicht der eng gefasste Begriff der Spek- 
tralverteilung einer von der Rontgenquelle der Apparatur emittierten Rontgenstrahlung zu verstehen, sondern ein erwei- 
terter BegrifF, der die unterschiedliche Gewichtung unterschiedlicher Spektralbereiche des Emissionsspektrums der 
RGntgenrohre auf Seiten der Rontgendetektoren beriicksichtigt. Ein gemessener Schwachungswert ergibt sich daher aus 25 
der direkten Schwachung des von der Rontgenrdhre emittierten Strahlenspektrums und dem spektralen Wirkungsgrad 
des verwendeten RSntgendetektors. Beide Werte sind anlagenspezifische GroBen und mussen entweder direkt oder indi- 
rekt mittels der Schwachungswerte von Eichproben ermittelt werden. Sie sind die Grundlage zur Berechnung der Isoab- 
sorptionslinien. 

[0043] In der Fig. 3a sind drei Verfahren 300 zur Modellierung bzw. zur Berechnung einer Schar von Isoabsorptions- 30 
linien skizziert, namlich eine theoretische Modellierung, eine experimentelle Bestimmung und eine theoretische Model- 
lierung mit einer Kalibrierung der Kurven durch experimentell bestimmte Parameter. 

[0044] Prinzipiell sind so viele Isoabsorptionslinien zu bestimmen, wie Schwachungswerte zum Abdecken der Spanne 
von Rontgenschwachungen in den Rontgenaufhahmen erforderlich sind. Dabei ist nicht fiir jeden theoretisch auftreten- 
den Schwachungswert eine Isoabsorptionslinie zu berechnen; nicht errechnete Isoabsorptionslinien konnen bei Bedarf 35 
durch Interpolation oder andere geeignete Mi ttelungs verfahren verfugbar gemacht werden. 

[0045] Die Grundschritte der theoretischen Modellierung sind im linken Ast des Ablaufschemas der Fig. 3a dargestellt. 
Im Schritt S302 werden zunachst die Daten der fiir eine Anlage spezifischen Rontgenemissionsspektren S(E) mit den 
yerfugbaren Rohrenspannungen als Parameter eingelesen. Die Spektralverteilungen der Rontgenstrahlung kdnnen hierzu 
im Vorfeld experimentell fiir jede einzelne Rontgenanlage ausgemessen werden, oder es werden die fiir einen speziellen 40 
Rontgenqueilentyp charakteristischen Daten verwendet Das Ermitteln der Detektorapparatefunktion w(E) erfolgt in 
Schritt S303. Auch hierzu kann im Vorfeld eine genaue Vermessung der Detektoranordnung vorgenommen werden oder 
aber es werden den Detektortyp charakterisierenden Daten wie z. B. dessen spektrale technische Spezifikation verwen- 
det. Die Berechnung der Isoabsorptionslinien in Form von Kurvenscharen Q(p, Z) bzw. Ui(p, Z) wird auf der Basis eines 
physikalischen Modells in Schritt S304 vorgenommen, das fur jede relevante {Combination von S(E) und w(E) die Ront- 45 
genschwachungen Q bzw. Ui fur Materialien mit unterschiedlichen Ordnungszahlen und bei unterschiedlichen Material- 
dichten nachbildet. 

[0046] Alternativ zur theoretischen Modellierung der Schritte S302 bis S304 konnen die Kurvenscharen der Isoabsorp- 
tionslinien auch experimentell ermittelt werden. Hierzu werden in Schritt S305 die Rontgenschwachungen von Eichma- 
terialien mit unterschiedlicher Dichte und mittlerer Ordnungszahl in der Rontgenapparatur bei verschiedenen relevanten 50 
Kombinationen von S(E) und w(E) gemessen. Die Messwerte bilden die Stiitzpunkte fur die folgende Berechnung der 
Kurvenscharen von Isoabsorptionslinien Q bzw. ui in Schritt S306. 

[0047] Als weitere Alternative konnen die auf theoretischer Basis modellierten Kurvenscharen Q bzw. Ui mit experi- 
mentell ermittelten Rontgenschwachungswerten kalibriert werden. In Schritt S307 werden die zum Eichen der theoreti- 
schen Kurvenscharen notwendigen Schwachungswerte wie oben fur Schritt S305 beschrieben mit geeigneten Eichmate- 55 
rialien bzw. Phantomen in der Rontgenanlage gemessen. Im Unterschied zur rein theoretischen Modellierung der Schritte 
S302 bis S304 ist bei diesem Verfahren die exakte Kenntnis der Rontgenemissionsspektren S(E) und w(E) nicht \foraus- 
setzung sondem Parameter der theoretischen Modellierung der Kurvenscharen von Isoabsorptionslinien Q bzw. Ui in 
Schritt S308. Das Kalibrieren der Kurven in Schritt S309 mit den in Schritt S307 experimentell ermittelten Eichwerten 
definiert schlieBlich Werte fur diese Parameter, die spezifisch fur die Rontgenemissionsspektren und Detektorapparate- 60 
funktionen der Rontgenapparatur sind. 

[0048] Mit der Ermittlung der Isoabsorptionslinien fiir die erforderlichen Rontgenschwachungswerte und Kombinatio- 
nen von S(E) und w(E) sind die Voraussetzungen fur eine Transformation von Bilddaten, die Schwachungswerte der 
Rontgenstrahlung beim Durchgang durch ein Gewebe reprasentieren in Bilddaten, die eine Verteilung der Ordnungszahl 
bzw. der Materialdichte im entsprechenden Gewebe reprasentieren geschafFen. 65 
[0049] Je nach Aufgabenstellung konnen die drei Verfahren zur Isoabsorptionslinienbestimmung auch gemischt ver- 
wendet werden. Beispielsweise konnen Werte, die experimentell nur ungenau oder nur mit groBem Aufwand oder gar 
nicht zu ermitteln sind, mithilfe einer theoretischen Modellierung erganzt oder in ihrer Genauigkeit prazisiert werden. 
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Die mit unterschiedlichen Methoden erschlossenen Daten werden dann in Schritt S3 10 zu einem einheitlichen Datensatz 
zusammengefasst und in Schritt S3 11 fiir die Bildtransformationen bereitgehalten. 

[0050] In der Fig. 3b ist das erfindungsgemafie Transformationsverfahren 320 dargestellt Es stiitzt sich auf die nach ei- 
nem der zuvor beschriebenen Verfahren 300 ermittelten und als Datensatz in Schritt S321 bereitgehaltenen Kurvenscha- 
5 ren von Isoabsorptionslinien. 

[0051] Eine Transformation erfolgt bildelementweise. Im folgenden wird von einer Transformation einer Rontgen- 
schwachungswertverteilung basierend auf zwei bei unterschiedlichen Rontgenstrahlenergiespektren aber identischer 
Aufnahmegeometrie aufgenommen RQntgenbildern ausgegangen. Dies ist die minimale \foraussetzung fur eine Durch- 
fuhrung einer erfindungsgemafien Iransformation. Jedoch konnen auch mehr als zwei Rontgenaufhahmen bei mehr als 

10 zwei unterschiedlichen Energieverteilungen der Rontgenstrahlung Verwendung finden. 

[0052] Die Auswahl eines zu transformierenden Bildelements wird im Schritt S322 getroffen und im folgenden Schritt 
S323 werden die Schwachungswerte Ci bzw. ui fiir dieses Bildelement aus dem ersten und C2 bzw. U2 aus dem zweiten 
Rontgenbild gelesen. Im anschlieSenden Schritt S324 erfolgt die Abfrage des fur die erste Rontgenaufhahme verwende- 
ten Rontgenstrahlspektrums Si(E) und der Detektorapparatefunktionen w L (E) sowie der entsprechenden Werte S 2 (E) und 

15 W2(E) fiir das zweite Rontgenbild. Diese Werte bilden die Parameter fur eine nachfolgende Auswahl der den jeweiligen 
Schwachungswerten zuzuordnenden Isoabsorptionslinien. Die Spektralverteilungen Si(E) bzw. Wi(E) konnen hierbei 
auch indirekt, z. B. iiber eine Abfrage der verwendeten Rohrenspannungen Ul bzw. U2 bzw. der Betriebsparameter der 
Rontgendetektoren ermittelt werden. 

[0053] Im Schritt S325 werden aus dem in Schritt S321 bereitgehaltenen Datensatz von Isoabsorptionslinien eine erste 
20 Kurve, welche die Bedingungen Ci bzw. Ux bei den Parametern Si(E) und Wi(E) erfullt und eine zweite Kurve, welche 
die Bedingungen C2 bzw. u 2 bei den Parametern S2CE) und w 2 (E) erfullt ausgewahlt. Ein Beispiel einer dergestalt erhal- 
tenen ersten 11 und zweiten 41 Isoabsorptionslinie ist im Diagramm 40 der Fig. 4 dargestellt. 

[0054] Der Schnittpunkt 42 als Schnittmenge beider Kurven 11 und 41 wird im Schritt S326 des Transformationsver- 
f aniens 320 berechnet. Der Kurvenschnitt 42 lasst sich z. B. durch eine lokale lineare Transformation oder mitteis itera- 
25 tiver Schnittpunktfindung ermitteln. Da die beiden Kurven 11 und 41 zwei unterschiedliche Schwachungswerte fur das 
selbe Bildelement und daher fiir einen identischen Teilbereich eines untersuchten Gewebes reprasentieren, mussen beide 
Schwachungswerte von der selben Material- bzw. Gewebeart verursacht sein. Die Koordinaten (p, Z) des Kurvenschnitt- 
punkts 42 geben daher die Materialdichte und die Ordnungszahl des dem Bildelement zuzuordnenden Gewebeteilbe- 
reichs wieder. 

30 [0055] SchlieBlich wird in Schritt S327 der so ermittelte Ordnungszahlwert Z in die Ordnungszahlverteilung als ent- 
sprechender Bildelementwert geschrieben, in Schritt S328 analog der ermittelte Materialdichtewert p in die Dichtever- 
teilung. Die Schritte S322 bis S328 werden fur alle verbleibenden Bildpunkte wiederholt, bis eine abschlieBende Bild- 
ausgabe in Schritt S329 erfolgen kann. Dabei kann der Schritt S324 iibersprungen werden, da die Spektralverteilungen 
Si(E) bzw. Wi(E) fur alle Bildelemente eines Bildes identisch sind. 

35 [0056] Die Genauigkeit der Bestimmung der Dichte p und der Ordnungszahl Z wird im wesentlichen vom Rauschan- 
teil der Schwachungswerte und der spektralen Auflosung der Messung beeinflusst Der Rauschanteil kann z. B. uber die 
verwendete Rontgendosis beeinflusst werden. Die spektrale Auflosung ist eine Funktion des Abstands und der Breite der 
zur Messung verwendeten Rontgenspektren. 

[0057] Generell erhalt man bei Rontgenspektren relativ niedriger Energie eine Bevorzugung der Rontgenschwachung 
40 durch den Photoeffekt, bei Rontgenspektren relativ hoherer Energie eine Bevorzugung der Rontgenschwachung durch . 
den Comptoneffekt. Genauer ausgedruckt ist der Einfluss der Ordnungszahl auf die Rontgenschwachungswerte einer 
Aufnahme bei geringerer Rontgenenergie relativ groBer als bei hdherer Rontgenenergie. Der Einfluss einer Material- 
bzw. Gewebedichte auf die Rontgenschwachungswerte verhalt sich dagegen genau umgekehrt. \brteilhaft wird daher zu- 
nachst ein erstes Rdntgenspektrum so gewahlt, dass ein deutlicher Anteil an den ersten Rontgenschwachungswerten vom 
45 Einfluss der Ordnungszahlen des untersuchten Gewebes oder Materials herruhrt und ein zweites Rdntgenspektrum wird 
dann so gewahlt, dass die Dichten des Untersuchungsobjekts einen deutlichen Einfluss auf die zweiten Rontgenschwa- 
chungswerte nehmen. 

[0058] Fiir die Computertomographie (CT) werden daher Quantenenergien der Rontgenstrahlspektren im Bereich zwi- 
schen 50 keV und 120 keV bevorzugt, so dass ein ausreichender Energieabstand zwischen einem ersten und einem zwei- 
50 ten Rdntgenspektrum gewahlt werden kann ohne die Rontgendosis in fur Patienten schadliche Bereiche erhohen zu miis- 
sen. In der Mammographie und der Mikro-CT konnen der untere Energiebereich sogar auf 10 keV bzw. 1 keV festgelegt 
werden. 

[0059] Die Material und Sicherheitspriifung verwendet i. A. hohere Rdntgendosen als bei der Untersuchung von Lebe- 
wesen vertretbar sind. Mit diesen hoheren Rontgendosen erhalt man ein meist wesentlich besseres Signal- zu-Rauschver- 

55 haltnis. AuBerdem konnen auch wesentlich hohere Rontgenstrahlenergien von bis zu 1 MeV verwendet werden. Damit 
kann, auch wenn die Ordnungszahl bei einer bestimmten Rontgenstrahlenergie nur zu einem geringen Anteil zur Ront- 
genabsorption beitragt, ihr Wert dennoch mit guter Genauigkeit bestimmt werden, sofern bei der verwendeten Rontgen- 
dosis der Rauschanteil der Rontgenabsorptionswerte deutlich geringer ist als der relative Anteil der Ordnungszahl an den 
entsprechenden Rontgenabsorptionswerten. Es kommt daher nicht so sehr auf die absolute energetische Lage der ver- 

60 wendeten Rontgenspektren an, sondern auf deren unterschiedliche energetische Lage zueinander. Die Wahl der absoluten 
Lagen wird dann von den verwendbaren Rontgendosen und den nachzuweisenden Stoffen beeinflusst. 
[0060] Ausreichende Ordnungszahl-Kontraste fur einen zufriedenstellenden Einblick in die chemische Zusammenset- 
zung eines Materials oder die biochemische Zusammensetzung eines Gewebes werden iiber die Rontgendosis und die 
verwendeten Rontgenspektren eingestellt. Die notwendige Auflosung in der Darstellung der Ordnungszahlverteilung 

65 hangt von der diagnostischen Aufgabenstellung ab. Fur eine quantitative Messung von typischen Korperbestandteilen 
wie z. B. Jod, Magnesium, Kohlenstoff, Stickstoff, Sauerstoff, Natrium, Phosphor, Schwefel, Chlor, Kalium, Mangan, 
Eisen, Kobalt, Kupfer, Zink, Selen, Molybdan oder dergleichen ist eine andere Z-Auflosung erforderlich als zum Erzie- 
len eines chemischen Kontrasts basierend auf z. B. dem Eisengehalt in einer Leber. Relativ gering sind die Anforderun- 
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gen an eine Z- Auflosung fiir ein Heraussegmentieren von Objekten aus einer Umgebung mit stark unterschiedlicher Ord- 
nungszahl. Prinzipiell muss die Z- Auflosung besser sein als die Differenz der noch zu unterscheidenden Ordnungszah- 
len. Gegenwartig ist eine Z- Auflosung von bis zu 0,1 erreichbar. 

[0061] Die Dichteverteilung kann zur Messung der Gewichtsverteilung eines Objekts gemaB der Gleichung 

m(r) = p(r) x Av(f ) (4) 5 

ausgewertet und zur Schwerpunktsbestimmung verwendet werden. Anwendungen sind z. B. die Dichtevermessung von 
Knochen oder des Gehirns, aber auch die Bestirnrnung des tatsachlichen Schwerpunkts eines Priifgegenstands. 
[0062] Zur Klassifizierung und Segmentierung von Korpermaterialien wie beispielsweise Leber, Nieren, weiBer und 10 
grauer Gehirnmasse etc. oder Materialklassen kann eine kombinierte Auswertung der Qrdnungszahl- und Dichtevertei- 
lung vorgenommen werden, indem eine Fensterung der Dichte und der iiber die Ordnungszahl bestimmten Gewebe- bzw. 
Materialart vorgenommen wird. Zusatzlich zum Segmentieren eines z. B. malignen Gewebes kann mit dieser Methode 
beispielsweise gleichzeitig dessen Stickstoffgehalt bestimmt und fiir die Dosisplanung in einer nachfolgenden Strahlen- 
therapie verwendet werden. 15 
[0063] Um einen Vergleich mit der klassischen, Schwachungswerte darstellenden Radiographie zu ermoglichen, kon- 
nen die urspriinglichen Bilder auch weiterhin als Standard-Rontgenbilder dargestellt werden. Falls diese nicht mehr vor- 
handen sind, lasst sich aus der Ordnungszahl- und Dichteverteilung eines Untersuchungsgegenstands riickwarts ein ab- 
wartskompatibles Standard Rontgenbild errechnen. 

[0064] Bei dem im Vergleich zur Computertomographie hohen Detektorrauschen von Rdntgenbildsystemen ist die Dy- 20 
nanuk der Rontgenbildkontraste und somit die erzielbare Z- und p-Auflosung reduziert. Unter Umstanden wiirde die 
Vorgabe einer hohen Z- bzw. p- Auflosung bei einer groBflachige Untersuchung eine unzulassig hohe Dosisbelastung be- 
deuten. Nach den "Obersichtsmessungen wird daher entsprechend den diagnostischen Erfordernissen in einem relativ 
kleinen Untersuchungsgebiet eine zweite Messung mit erhohter Dosis vorgenommen. In der Mammographie empfiehlt 
sich diese Methode zur Untersuchung von Mikro-Calcifizierungen. 25 
[0065] In der Fig. 5 sind drei verschiedene ROntgensysteme zur Illustration dreier Grundprinzipien von Rontgendurch- 
leuchtungsverfahren dargestellt. 

[0066] Fig, 5a zeigt ein Rontgenbildsystem 500 zur projektiven Abbildung eines Untersuchungsobjekts 503 auf einen 
zweidimensionalen Flachendetektor 502. Ein Patient 503 oder ein Gegenstand 503 befinden sich darin im allgemeinen 
wahrend einer Rontgenabbildung relativ zur Anordnung aus Rontgenrohre 501 und Detektor 502 in Ruhestellung. 30 
[0067] Die Fig. 5b gibt als Beispiel einer Computertomographie- Vorrichtung schematisch das Prinzip eines Facher- 
strahlgerats 500 mit rotierendem Detektorsystem 502 wieder. Die von der Rontgenrohre 501 facherformig emittierten 
Rontgenstrahlen durchleuchten das Untersuchungsobjekt 503 aus einer bestimmten Winkelstellung und treffen schlieB- 
lich auf eine zeilenformige Anordnung diskreter Detektoren. Ein Aufhahmezyklus umfasst eine Vielzahl solcher Durch- 
leuchtungen in unterschiedlichen Winkelstellungen zum Untersuchungsgegenstand. 35 
[0068] Bei dem in der Fig. 5c illustrierten Typ einer Rontgendurchleuchtungsvorrichtung 500 wird der Untersuchungs- 
gegenstand 503 mit einer TYansportvorrichtung wie z. B. einem TVansportband 506 durch den Strahlengang 504 der An- 
lage gefuhrt. Dieser Typ von Rontgenapparatur wird hauptsachlich in der Sicherheitspriifung verwandL Als Detektor 
wird ein linienformiger, eindimensionaler, oft L-formiger Detektor verwendet. Beim Durchfahren durch den Strahlen- 
gang 504 wird der Gegenstand 503 von den von der Rontgenrohre 501 emittierten Rontgenstrahlen zeilenweise abgeta- 40 
stet und aus den so erstellten Bildzeilen ein vollstandiges Durchleuchtungsbild des Gegenstands zusammengesetzt 
[0069] Die Fig. 6 zeigt verschiedene Moglichkeiten der Beeinflussung des efFektiven Rontgenstrahlspektrums, wie es 
sich aus dem Zusammenwirken der von der Rontgenrohre 501 emittierten Strahlungsverteilung 610 (611) und der spek- 
tralen Charakteristik der Detektorapparatefunktion 612 (620) ergibt. 

[0070] Fig. 6a illustriert ein schnelles Umschalten der Rohrenspannung von einem Wert Ui auf einen Wert U2, so dass 45 
ein Wechsel zwischen zwei Rontgenemissionsspektren 610 und 611 erzielt wird. Vorteilhaft wird jede Teilaufhahme ei- 
nes Aufhahmezyklus zuerst bei einer der beiden Rohrenspannungen und unmittelbar anschliefiend bei der anderen der 
beiden Rohrenspannungen aufgezeichnet um Artefakte durch Bewegungen z. B. eines Patienten wahrend des Aufnah- 
mezyklus zu rninimieren. Bei Rontgenapparaten 500 mit kurzen Messzeiten wie z. B. Rontgenprojektionssystemen kann 
auch jeder Aufnahmezyklus vollstandig bei jeweils einer der beiden Rohrenspannungen durchfahren werden. Statt einem 50 
Umschalten der Rohrenspannung konnen auch Filter oder Filtersysteme zwischen Rontgenrohre und Untersuchungsob- 
jekt bzw. zwischen Untersuchungsobjekt und Rontgendetektor zur Beeinflussung des Rontgenspektrums eingebracht 
werden. Der dargestellt Betrieb der Rontgenrohre bei einer Betriebsspannung Ui wird dann ersetzt durch einen Betrieb 
bei einer Standard-Betriebsspannung und Verwendung von keinem oder einem ersten Filter bzw. Filtersystem; ein Be- 
trieb bei U2 wird durch Verwenden eines zweiten Filters bzw. Filtersystems ersetzt. 55 
[0071] Da eine Rontgenrohre 501 keine monoenergetische Strahlung sondern ein relativ breites Spektrum von Ront- 
genstrahlung 504 emittiert, kann man, wie in der Fig. 6b gezeigt in einem Messvorgang mehrere Aufnahmen bei unter- 
schiedlichen Rontgenspektren erhalten indem man einen energiesensitiven Detektor 620 verwendet. Dieser liefert fur 
i. a. nebeneinander angeordnete Spektralbereiche getrennte Messsignale. Man erhalt Schwachungswerte somit simultan 
fiir unterschiedliche, voneinander getrennte Spektralbereiche des Rontgenspektrums also in einem Aufhahmezyklus 60 
eine, von der Ausfuhrungsform und Beschaltung des Detektors definierte Anzahl von Rontgenbildem bei unterschiedli- 
chen Strahlenergien. Solche Detektoren 620 konnen als Schichtaufbau-Detektoren realisiert werden, wobei ausgenutzt 
wird, dass die Eindringtiefe von Rontgenstrahlung in das Schichtsystem des Detektors von der Energie der Rontgenquan- 
ten bestimmt wird. Altemativ zu Schichtdetektoren konnen Quantenzahler als energiesensitive Detektoren verwendet 
werden. 65 
[0072] Eine dritte Variante zur Beeinflussung des Rontgenemissionsspektrums zeigt die Fig. 6c. Hierbei wird ein Ma- 
terial 630 eines aufgabenspezifischen Elements z. B. in Form schnell drehender Filterrader 630 zwischen Rontgenrohre 
501 und Untersuchungsobjekt 503 oder zwischen Untersuchungsobjekt 503 und Rontgendetektor 502 platziert. Jedes 
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Element absorbiert auf eharakteristische Weise Rontgenstrahlung in einem ebenfalls charakteristischen Spektralbereich. 
Befindet sich ein Fiugel des Filterrads im Strahlengang, so wird ein fur dieses Element typischer Spektralbereich im 
Strahlengang gedampft bzw. im giinstigsten Fall eliminiert. Fiir jede Teilaufiiahme eines Aufhahmezyklus werden nun 
eine Messung ohne und eine Messung mit einem FilterradflUgel im Strahlengang vorgenommen. Ein dem Filterradele- 
ment entsprechendes Element im Untersuchungsobjekt kann dann nur in den Messungen ohne Fllterradabdeckung des 
Strahlengangs einen Rontgenschwachungswert 632 aufweisen. Im anderen Fall ist der gemessene Schwachungswert 633 
idealerweise gleich Null. Auf diese Weise kann mit hoher Prazision das Vorhanden- bzw. Nichtvorhandensein eines be- 
stimmten Elements, wie beispielsweise Eisen oder Kupfer in einem Untersuchungsobjekt nachgewiesen werden. Diese 
Nachweismethode ist insbesondere in der Mammographie und der Materialprufung von Bedeutung. 

Patentanspriiche 

1. Verfahren zur Ermittlung der Verteilungen von Dichte und Ordnungszahl in einem Untersuchungsobjekt (503) 
mit Schritten 

zum Aufzeichnen einer ersten Verteilung von Rontgenabsorptionswerten des Untersuchungsobjekts (503) bei ei- 
nem ersten Rdntgenspektrum und zum Aufzeichnen zumindest einer zweiten Verteilung von Rontgenabsorptions- 
werten des Untersuchungsobjekts (503) bei einem zweiten Rontgenstrahlspektrum, 

zum Bestimmen einer ersten funktionalen Abhangigkeit (11) eines ersten Rontgenabsorptions wertes der ersten Ver- 
teilung von Rontgenabsorptionswerten von Dichte und Ordnungszahl und 

zum Bestimmen zumindest einer zweiten funktionalen Abhangigkeit (41) eines, dem ersten Rontgenabsorptions- 
wert zugeordneten zweiten Rontgenabsorptionswertes der zweiten Verteilung von Rontgenabsorptionswerten von 
Dichte und Ordnungszahl, und 

zum Ermitteln eines Werts fur die Dichte und die Ordnungszahl aus einem Vergleich der ersten funktionalen Abhan- 
gigkeit (11) mit der zweiten funktionalen Abhangigkeit (41) und/oder weiteren funktionalen Abhangigkeiten. 

2. Verfahren nach Anspruch 1, dadurch gekennzeichnet dass die Bestimmung der funktionalen Abhangigkeit (11, 
41) der Rontgenabsorptionswerte von Dichte und Ordnungszahl fur zumindest ein Rontgenstrahlspektrum mittels 
Referenzmessung an einer Eichprobe erfolgt. 

3. Verfahren nach Anspruch 1 oder 2, dadurch gekennzeichnet dass die Bestimmung der funktionellen Abhangig- 
keit (11, 41) der Rontgenabsorptionswerte von Dichte und Ordnungszahl fur zumindest ein Rontgenstrahlspektrum 
in Form einer Simulation auf der Basis eines physikalischen Modells erfolgt. 

4. Verfahren nach Anspruch 1, 2 oder 3, dadurch gekennzeichnet dass ein Umformen der \ferteilungen der Ront- 
genabsorptionswerte in eine Verteilung der Dichte und eine Verteilung der Ordnungszahl fur jeden der zugeordne- 
ten Rontgenabsorptionswerte der ersten Verteilung der Rontgenabsorptionswerte und der weiteren Verteilungen der 
Rontgenabsorptionswerte auf der Grundlage der Ermittlung eines Wertepaares (42) fur Dichte und Ordnungszahl so 
vorgenommen wird, dass das Wertepaar (42) die bestimmten funktionalen Abhangigkeiten der Rontgenabsorption 
von Dichte und Ordnungszahl fiir das erste Rontgenstrahlspektrum und zumindest ein wei teres Rontgenstrahlspek- 
trumerfiillt. 

5. Verfahren nach einem der Anspruche 1 bis 4, dadurch gekennzeichnet dass das erste Rdntgenspektrum eine 
Quantenenergie aufweist, die relativ zur Quantenenergie des zweiten Rontgenspektrums eine Rontgenabsorption 
durch den Photoeffekt begunstigt. 

6. Verfahren nach einem der Anspruche 1 bis 5, dadurch gekennzeichnet dass zum Verandern eines Rfintgenstrahl- 
spektmms fiir das Aufzeichnen des Untersuchungsobjekts (503) eine Veranderung zumindest eines Betriebspara- 
meters der Rontgenrohre und/oder eine Veranderung der Detektorcharakteristik vorgenommen wird. 

7. Verfahren nach einem der Anspruche 1 bis 6, dadurch gekennzeichnet dass das Aufzeichnen zumindest einer 
Teilregion des Untersuchungsobjekts (503) mit erhohter Rontgendosis erfolgt. 

8. Rontgenapparatur (500) zur Ermittlung der Verteilungen von Dichte und Ordnungszahl in einem Untersuchungs- 
objekt mit 

einer Rontgenquelle (501) zur Emission von Rontgenstrahlung (504), 

einem Rontgendetektor (502) zum Nachweis der von der Rontgenquelle emittierten Rontgenstrahlung (504) und 
zum Umsetzen der Rontgenstrahlung (504) in elektrische Signale fur eine weitere Verarbeitung, und 
einer Signalverarbeitungseinrichtung (505) fur das Verarbeiten der elektrischen Signale des Rontgendetektors (502) 
dadurch gekennzeichnet, 

dass die Signalverarbeitungseinrichtung (505) die Verteilungen von Dichte und Ordnungszahl in einem Untersu- 
chungsobjekt (503) auf Basis zumindest zweier, an der Rontgenapparatur (500) mit unterschiedlichen Rontgen- 
strahlspektren aufgezeichneten Verteilungen einer Rontgenabsorption des Untersuchungsobjekts (503) nach einem 
der Verfahren der Anspruche 1 bis 7 ermittelt. 

9. Rontgenapparatur nach Anspruch 8, dadurch gekennzeichnet, dass die Rontgenquelle (501) in einem ersten Be- 
triebszustand ein erstes Rontgenstrahlspektrum (610) und in einem zweiten Betriebszustand ein davon verschiede- 
nes zweites Rontgenstrahlspektrum (611) emittiert. 

10. Rontgenapparatur nach Anspruch 8 oder 9, dadurch gekennzeichnet, dass der Rontgendetektor (502) spektrale 
Teilbereiche der von der Rontgenquelle (501) empfangenen Rontgenstrahlung (504) in voneinander unabhangige 
elektrische Signale umsetzt 

11. Rontgenapparatur nach Anspruch 8, 9 oder 10, dadurch gekennzeichnet, dass zum Aufzeichnen einer ersten 
Verteilung einer Rontgenabsorption des Untersuchungsobjekts (503) ein Material (630) in den Strahlengang der 
Rontgenstrahlung (504) eingebracht wird und zum Aufzeichnen einer zweiten Verteilung einer Rontgenabsorption 
des Untersuchungsobjekts (503) kein entsprechendes Material in den Strahlengang der Rontgenstrahlung (504) ein- 
gebracht wird. 

12. Computersoftwareprodukt zur Ermittlung der Verteilungen von Dichte und Ordnungszahl in einem Untersu- 
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chungsobjekt (503) auf Basis von zumindest zwei bei unterschiedlichen Rontgenstrahlspektren aufgezeichneten 
Bilddatensatzen der Verteilung von Rontgenab sorptions werten in dem Untersuchungsobjekt (503) nach einera der 
Verfahren 1 bis 7, 
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